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ABSTRACT 

Adams et al. (1997b) have noted that according to our current understanding of the 
unification of fundamental interactions, there should have been phase transitions asso- 
ciated with spontaneous symmetry breaking during the inflationary era. This may have 
resulted in the breaking of scale-invariance of the primordial density perturbation for 
brief periods. A possible such feature was identified in the power spectrum of galaxy 
clustering in the APM survey at the scale k ~ 0.1/iMpc -1 and it was shown that 
the secondary acoustic peaks in the power spectrum of the CMB anisotropy should 
consequently be suppressed. We demonstrate that this prediction is confirmed by the 
recent Boomerang and Maxima observations, which favour a step-like spectral feature 
in the range k ~ (0.06 — 0.6)/iMpc -1 , independently of the similar previous indication 
from the APM data. Such a spectral break enables an excellent fit to both APM and 
CMB data with a baryon density consistent with the BBN value. It also allows the 
possibility of a matter-dominated universe with zero cosmological constant, which we 
show can now account for even the evolution of the abundance of rich clusters. 



1 INTRODUCTION 



It is commonly assumed by astronomers that inflation pre- 
dicts a scale-invariant 'Harrison-Zeldovich' (H-Z) spectrum 
of scalar density perturbations, S 2 1 (k) oc fc n_1 with n = 1; 
for example this is a standard input in calculations of the 
expected large-scale structure (LSS) and cosmic microwave 
background (CMB) anisotropy. In fact in the usual 'slow- 
roll' inflationary scenario there is a gradual steepening of the 
spectrum with increasin g k (decreasing scale), as the end of 



(sCDM) model with LSS data, in particular the power spec- 
trum of galaxy clustering, P{k), the abundance of rich clus- 
ters quantified by the parameter as, and large-scale stream- 
ing velocities, <r„(40 /i _1 Mpc). We emphasise that this model 
is however ruled out by the same data if a scale-inv ariant 
spectrum is assumed (Efstathiou, Bond & White 1992). Such 



conclusions are clearly not robust given that we have as yet 
no standard model of inflation. 



inflatioiriri approached (Mrrkhanov fe ChlblHov r981r| ff. 
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Subsequently it was realized (Adams, Ross & Sarkar 



1997b) that the spectrum may not be even scale-free because 



hardt & Turner 1983). This cannot be adequately modelled 
by a spectrum with constant "tilt" (n < 1) as is occasion- 
ally adopted, since the steepening is scale- dependent for any 
polynomial inflaton potential.^ Even such small departures 
from scale-invariance can be quite significant, for example in 
a supergravity model with a cubic inflaton potential, Adams, 
Ross & Sarkar (1997a) found n ~ (N — 2)/(N + 2), where 
N ~ 51 + In(fc -1 / 3000 h^Mpc) is the number of e-folds 
of expansion from the end of inflation. This gradual sup- 
pression of small-scale power was shown to be adequate to 
reconcile the COBE-normalised standard cold dark matter 



* The index n is scale-independent only for an exponential po- 
tential ('power-law inflation'); moreover n can be very close to 
unity if inflation ends not through the steepening of the inflaton 
potential but, for example, due to the dynamics of a second scalar 
field ('h ybrid inflation''!. The scale-dependence of n in various in- 



the rapid cooling of the universe during primordial inflation 
can result in spontaneous symmetry breaking phase transi- 
tions which may interrupt the slow roll of the inflaton field 
for brief periods. This is in fact inevitable in models based on 
N — 1 supergravity, the phenomenologically successful ex- 
tension of the Standard Model of particle physics and the ef- 
fective field theory below the Planck scale. During inflation, 
the large vacuum energy breaks supersymmetry giving oth- 
erwise massless fiel ds ('flat directions') a mass of order the 



Hubble parameter (Dine et al. 1984; Soughlan et al. 1984), 
causing them to evolve rapidly to the asymmetric global 
minima of their scalar potential. Such 'intermediate scale' 
fields are generic in models derived from superstring/M- 
theory and have gauge and Yukawa couplings to the thermal 
plasma so are initially confined at the symmetric maxima of 
their scalar potentials. Consequently it takes a (calculable) 
finite amount of cooling before the thermal ba rrier disap- 
pears and they are free to evolve to their minima ( Yamamotc 



national y models has been reviewed by Lyth & Riotto (1999). 



1986 ; Barreiro et al. 1996 ) . When a symmetry breaking tran- 
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sition occurs, the mass of the infiaton field changes suddenly 
(through couplings in the Kahler potential), temporarily vi- 
olating the slow-roll conditions and interrupting inflation 
(Adams et al. 1997b) .£] Thus the density perturbation is ex- 
pected to have a (near) H-Z spectrum for the first ~ 10 
c-folds of expansion followed by one or more sudden depar- 
tures from scale-invariance lasting ~ 1 e-fold. In order for 
such spectral features to be observable in the LSS or CMB, 
it is of course necessary that they occur within the last ~ 50 
e-folds of inflation, corresponding to spatial scales going up 
to the present Hubble radius Hg 1 ~ 3000 /i _1 Mpc. Since 
the density perturbation can be observed on scales from the 
Hubble radius down to ~ 1 ft Mpc, corresponding to about 
8 e-folds of expansion, it would be not unreasonable to ex- 
pect at least one such spectral break to be seen today. 

Motivated by this an attempt was made to recover the 
primordial perturbation spectrum from extant observations 
of LSS and CMB ani sotropy, specifically the APM survey 
( Maddox et al. 1990a ) and the COBE observations (Smoot 
et al. l f>92| ). We recall that the spectrum of rms mass fluc- 
tuations at the present epoch (per unit logarithmic interval 
of wavenumber fc) is 



below unity) at fc ~ 0.1 ftMpc . The sharp drop occurs over 
~ 1 e-fold of expansion as is indeed expected in the phase 
transition model. This feature in Papm(&) had been noted by 
Baugh & Efstathiou (1993) and was also identified by Pea- 
cock (1997) in the power spectrum of IRAS galaxies after 
correcting for redshift-space distortions. Given that the fea- 
ture appears on the scale of the Hubble radius at the epoch 
of (dark) matter domination, it is perhaps natural to inter- 
pret it as reflecting a departure from the sCDM paradigm 
as these authors did. However Adams et al. (1997b) pursued 
the alternative possibility that it is in fact present in the pri- 
mordial spectrum and asked what would be the expectation 
for the angular anisot ropy of the CMB i f this were so. Using 
the COSMICS code ( |Bertschinger 1995| ) they calculated its 
angular power spectrum and found that the height of the 
secondary acoustic peaks is suppressed by a factor of ~ 2 
relative to the prediction of the standard COBE-normalised 
sCDM model which assumes a scale-free H-Z spectrum. 

We reexamine this predic tion in the light of recen t ob- 
servations by the Boomerang (|dc Bcrnadis et al. 2o"oC ) and 
MAXIMA- 1 ( |Hanany et al. 2000| ) experiments which do in- 
dicate a significant suppression of t he second acoustic peak 



A 2 (k) k 



k J P(k) 



9-h 



It has been noted by ma n y authors ( Tegmark fc Zaklarriagc 
200C|; |Lange et al. 200l|; |Ba.lbi et al. 2000|; |Hu et al. 200C 
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Jaffe et al. 2000 



where t |he density perturbation is evaluated at the present 
Hubble radius, i.e. at k — Ho. The (dimensionless) matter 
' transfer function' for CDM models can be approximated by 
(iBond & Efstathiou 1984)) 



Esposito et al. 2001 



Padmanabhan & Sethi 



TcDM(fc) 



1 + {ak+ (bk) :i/2 + (cfc) 2 }' 



-i/» 



(2) 



with a = 6.4P" 1 /i _1 Mpc, b = 3r _1 ft _1 Mpc, c = 
1.7r _1 fe -1 Mpc and v — 1.13. H ere the 'shape parameter' is 



defined as (Bond fc Jaffe 1999) 



T = fl Jjg— [ n B(l+VaK/O m )— 0.08] 



(3) 



where h = 77o/100Kms _1 Mpc _1 is the Hubble parameter 
and fl m — ficDM + fiB, S^cdm (^b) being the fraction of the 
critical density in cold (baryonic) matter. Thus T = 0.5 for 
sCDM which assumes Qb = 0.03, Sl m = 1 and h — 0.5. 

Following Baugh & Gaztanaga (1996), the three- 
dimensional PAPM(fc) inferred from the angular correlation 
func tion of galaxies in the APM survey ( Baugh fc Efstathiou 
1993) was subj ected to an inversion p rocedure ( |jah*i7^c 



& Whi e 1995 Peacock fc Dodds 1996) to recover the lin- 



ear spectrum of fluctuations and obtain its spectral index, 
dlnP(fc)/dlnfc, as a function of fc. Assuming that any bias 
between APM galaxies and dark matter is scale-independent, 
the primordial spectral index n(fc) was then obtained simply 
by subtracting off the slope of the transfer function (^) for 
each value of fc (Adams et al. 1997b). 

The key finding was that nAPM(fc) shows a dramatic 
departure from the H-Z value of n = 1 in the range fc ~ 
(0.05 — 0.6)h Mpc -1 , dropping to a value around zero (~ 3a 



t When (re)heating occurs at the end of inflation such fields may 
again be forced back to the symmetric maximum, undergoing 
symmetry breaking a second time when the universe cools down to 
the electroweak scale in the radiati on-dominated era and driving 
a late phase of 'thermal inflation' ( |Lyth &: Stewart 1996[ ). 



200C ) that this requires a high baryon density, significantly 
above the value required by light element abundances fol- 
lowing from standard Big Bang nucleosynthesis (BBN). As 
the latter is widely considered to be one of the "pillars" of 
modern cosmology, it is important to ask whether this con- 
flict may not be resolved by relaxing the assumption made 
in these analyses that the primordial spectrum of perturba- 
tions is scale-free. 



2 RECONSTRUCTING THE PRIMORDIAL 
SPECTRUM 

We begin by describing the range of values that we consider 
for the cosmological parameters. This represents the set of 
priors, resulting from other observations, that we adopt for 
the reconstruction of the primordial spectrum. In §2.2, we 
discuss how LSS data can be used to recover the primordial 
spectrum and predict the CMB anisotropy (§2.3). In §2.4 we 
present a simple parameterization for the primordial spec- 
trum and use this to fit the new CMB as well as APM data. 
Finally in §2.5 we confront the observed abundance of rich 
clusters with the expectations in our model. 



2.1 Cosmological parameters 

Let us first briefly discuss the choice of cosmological param- 
eters. As mentioned before, F = 0.5 in the sCDM m odel, but 
r ~ 0.2 gives in fact a better fit to the APM date 



et al. 1990b; 


Efstathiou, Bond & White 1992| 


stathiou 199J 


; Peacock & Dodds 1994J) and c 



e.g. in a low density model with f2 m ~ 0.4.^] The alternative 



t However structure formation does not necessarily require a low 
density universe as is often claimed (e.g. Bahcall et al. 1999); for 
example with a gradually tilting primordial spectrum obtained 
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ibi lity of a Hubble parameter as low as h ~ 0.3 ([Bartlett the mass fraction Y( 4 He) =0.243 ± 0.003 inferred ( [izotov 



possih 

et al. IS 95) is definitively ruled out by the HST Key Project 



data on 25 galaxies which indicates h — 0,71±0.06 ( |Mould et 
al. 200(||h_ Weight ng such direct observations against other 



Thuan & Lipovetsky 1998 ) from observations of metal-poor 



fundan ental physics approaches, Primack (200(1) advocates 



h — 0.65 ±0.08. However a value as low as 0.45 ( jParodi et al, 
200C |) o| r as high as 0.9 (Ibnry et al. 2000) still seems quite 
possibl e so the quoted uncertainty cannot be interpreted as 



blue compact galaxies (BCGs), and the 'Spite pla teau' of 
7 Li/H=(1.73 ± 0.12 ) x 10~ 10 seen in Pop II stars flBonifa- 
cio fe Molaro 199%. Note however that an improved fit is 



obtained ( |Fiorcntini et al. 1998 ) at a significantly smaller 



auss ian standard deviation. Instead wc assume that all 



a ga 

values in the following "2a" range are equally probable: 



value of fl B = (0.0065l"o.ooM4) /l ~ : ' if one ad opts the 'high' 
value of D/H=(1.9 ± 0.4) x 10~ 4 claimed ( gongaila et al 
1994) in another QAS [ but disputed subsequently (Buries 



0.49-0.81. 



Kirkman fe Tytler 199S), although consistent with a possi- 



ble high value ( Tytler et al. 199E ) in yet another QA SJ, and 
( 4 ) the 'low' value of Y( 4 He)= 0.234 ± 0.005 obtained jOlive ' 



The current preference is in fact for a cosmological constant 



Skillman & Steigman 1997) from analysis of older BCG 



dominated A at universe gritb Qa — aTirl ° m — 1/3 ( p g 



data [and supported by ano ther recent observation ( Peim 



Bahcal ll et al 1999). The Boomerang and MAXTMA oh- 



servations of the fir st acoustic peak indicate that the uni- 
verse is indeed flat (Jaffe et al. 2000), but we will consider 



other ACDM models with values of 57a ranging down to 
zero, maintaining fim ± Qa = 1. This is motivated by our 
concern that the only direct evidence suggesting Qa ~ 0.7 
is the small curvature of the Hubble diagram for T y pe la 
supern ovae at redshift z ~ 1 ( Perlmutter et al. 199£ ; 



et al. 1 )9£ ) which may be subject to unrecognised system- 
atic effects. (We note that other arguments for a low Q, n 
fro m th e evolution of the cluster number density (like" et al 



1998) or the clustering in the Lyman-a forest (Weinberg et 
al. 199! ]) rest on the assumption of a scale-free spectrum.) 
Of course for any particular choice of (f2 m ,f2A), one has to 
check that the age of the universe excee ds the inferre d age 
of globular clusters, t G c = 12.8 ± 1 Gyr ( [Krauss 200o| ). 

Another important parameter for the present study is 
the BBN value of Ob- This has undergone substantial revi- 
sion since the first studies of sCDM so we briefly review the 
present situation. The baryon density which provides a good 
fit to the inferred pri mordial abundances o f D, 4 He and 7 Li 
in standard BBN is (Fiorcntini et al. 1998) 



(o.oi9tr »- 



(•») 



taking all (correlated) reaction rate uncertainties into ac- 
count.^] This fit is based on the Tow' value of D/H=(3.4 ± 
0.25) x 10~ 5 found ( [Buries fe Tytler 1998| ) in two quasar ab- 
sorption system s (QAS) at high redsh ift [and subsequently 
in a third one (Kirkman et al. 200C)], the 'high' value of 



from 'new' supergravity inflation, the small-scale power is suffi- 
ciently suppressed for a f2 m = 1 CDM model to be compatible 
with LSS data (Adams et al. 1997a, see also White et al. 1995). 
Such 'tilt' can preferentially suppress the second acoustic peak if 
the baryon fraction is high (Adams et al. 1997a, see also White 
et al. 1996). Several groups have studied the implications of tilt 
for the recent CMB data (e.g. Kinney, Melchiorri & Riotto 2000, 
Covi & Lyth 2000, Tegmark, Zaldariagga & Hamilton 2001). 
§ Subsequently Buries et al. (1999) suggested that some of the 
reaction rate uncertainties had been overestimated; thus they 
inferred a somewhat tighter range Qb = (0.019 ± 0.0009)/i — 2 . 
Izotov et al. (2000) found however a better fit with Qb = 
{0m&t^ ° \l)h~ 2 , taking Y( 4 He)=0.2452±0.0015 from observa- 
tions of the two most metal-deficient BCGs, D/H=(4.35±0.43) X 
10 — 5 from the analysis of O AS data using kinematic models incor - 
porating 'mesoturbulence' (Lcvshakov, Kegel & Takahara 199E), 
and 7 Li/H=(2.24±0.57) X 10~ lu allowing for some deple tion due 
to rotation in Pop II stars ( Vauclair fe Charbonnel 199ij ). 



bert, Peimbert fe Ruiz 2000)]. An independent lower limit of 



Qb > 0.0125ft (Weinberg etal. 1997) set by observations of 



the Lyman-a forest favours the higher value in Equation (^), 
which also allows a better fit to the observed distribution of 
line widths in comparison with model hydrodynamic simula- 
tions (e.g. Theun s et al. 1999). An audit of luminou s matter 
in the universe (Fukugita, Hogan & Peebles 1998) gives a 



much weaker lower limit of Qb > .007. A conservative up - 
per limit on the baryon density is (Kernan & Sarkar 1996) 



fi B < 0.033 /T 



(6) 



if the primordial deuterium abundance is bounded from be- 
l ow by its typical i nterstellar value D/H=(1.5 ± 0.2) x 10 -5 
(McCullough 1992). Therefore we will consider a wide range 
of possible values in our analysis, but adopt Ob = 0.045, 
h = 0.65 as the default values (so that ficDM = O m — 0.045) 
when not stated otherwise. 



2.2 Primordial spectrum from the APM data 

There are several arguments (reviewed in Appendix A) that 
on scales O.OlftMpc -1 < k < 0.6ftMpc _1 , which are at 
most wea kly non-linear, the APM galaxy power spectrum 
-Papm (fe) ( Baugh fe Efstathiou 1993 ) is an unbiased (or mod- 
erately linearly biased) tracer of the mass. The linear power 
spectrum recovered under this assumption from Papm(&), 
using the method of Jain et al. (199 5), is well fitted in this 
range by (Baugh & Gaztafiaga 1996): 



flta(fc) 



7 x 10 5 fc(ft- 1 Mpc) 

[1 + (fc/fcc) 2 ] 1 - 6 



(V 



where k c = 150(if /c) ~ 0.05 ftMpc -1 . We will use the 
common convention of expressing the matter power spec- 
trum as the primordial spectrum of matter fluctuations 
P°(k) oc k6ft(k), folded with the transfer function (|2j): 



Plin(fc) 



P°(k) T 2 {k) =Ak n T 2 (k), 



(8) 



where the (dimensionful) normalisation constant A is deter- 
mined by the large-scale CMB anisotropy. [F or sCDM with 
n = 1 normalised to the COBE 4-year data ( Bennett et al. 
1996; ), 5n = 1.92 x 10" 5 with 9% uncertainty ( Bunnfe White 
1997) so A = (5.9 ± 1.1) x 10 5 ( ft _1 Mpc) 4 .] 

Given some estimate of the (linear) matter power spec- 
trum, as in Equation ([?]), the standard approach is to as- 
sume a scale-invariant initial spectrum, i.e. P°(k) cx k, and 
extract the value of F that best fits the resulting constraint 
on T(k). Alternatively one can extract P°(k)(= P lin /T 2 ) 
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Figure 1. Reconstruction of the primordial density power spec- 
trum from the APM data, adopting a CDM shape parameter V of 
0.5 (dot-dashed line) and 0.2 (dashed line). A Harrison-Zeldovich 
spectrum (full line) is shown for comparison. 



for a specific choice of V, i.e. for a given transfer function. 
We apply this latter approach to the APM data by simply 
taking the ratio of Equation (^) to Equation (Q), obtaining: 



P°(k) 



A x k 
Ak 
A 2 k 



[l + ( a fc+(6fc) 3 / 2 + (cfc) 2 )"] 2/ " 

[i+(fc/MT 6 



k < ki, 

ki <k< ki, (9) 
fca < k, 



where A = 7 x 10 5 ( ft~ 1 Mpc) 3 and we have adopted a suit- 
ably normalised H-Z spectrum outside the range (k < ki = 
0.01 h Mpc , k > &2 = 0.6ftMpc _1 ) accessible to galaxy 
clustering observations. 

Figure |l| shows the recovered spectrum ^ for two 
choices of T corresponding to the sCDM model and a low 
density variant. Note that the F = 0.5 reconstruction has 
significantly less power than a scale-free H-Z spectrum on 
scales k > 0.1 h Mpc -1 , while the F = 0.2 reconstruction 
is closer to a H-Z spectrum but has relatively more power. 
However as we shall see the latter possibility does not give 
a good fit to the Boomerang/MAXIMA data with the BBN 
value (^) of the baryon density. To ensure compatibility with 
BBN requires a break in the primordial spectrum. 




200 400 600 
rnultipole 



800 



Figure 2. The predicted CMB angular power spectra correspond- 
ing to the reconstructed primordial spectra for T = 0.5 (dot- 
dashed line) and T = 0.2 (dashed line). The corresponding re- 
sults for a primordial H-Z spectrum is also shown for T = 0.5 
and r = 0.2 (thick and thin full lines). The data are from the 
Boomerang (crosses) and MAXIMA (diamonds) experiments. 



r = 0.2 case (corresponding e.g. to f2 m = 0.35, £Ia = 0.65), 
the inferred primordial spectrum is a factor of ~ 2 above 
the H-Z one at k ~ 0.04 ft Mpc -1 so the CMB anisotropy 
should be proportionally boosted at I > 250. 

To chec k these expectations, we ran the Boltzmann code 
CMBFAST flSeljak fc Zaldarriaga 1996| ) with the above re- 
constructed primordial spectra as well as a H-Z spectrum 
with the same cosmological parameters. As seen in Figure 
^, the angular power resulting from the APM inversion with 
r = 0.5 is indeed depressed below the H-Z case from the 
second acoustic peak onwards, as noted by Adams et al. 
(1997b). By contrast, for the T — 0.2 reconstruction all the 
acoustic peaks are boosted above the H-Z case. It is clear 
that the recent CMB observations favour a high density 
universe. Moreover the observed suppression of the second 
acoustic peak in the Boomerang/MAXIMA data confirms 
the prior expectation for a primordial density perturbation 
with broken scale-invariance. To quantify this we now per- 
form detailed fits to the data. 



2.3 Acoustic peaks 

In linear theory the co-efficents Ci in the spherical harmon- 
ics expansion of the CMB anisotropy are a linear projection 
of P{k), so the effect of changing P°(k) would just be to 
proportionally alter the angular power for the multipoles I 
corresponding to the relevant values of k. We can approxi- 
mate I ~ 7r/#LSS where #lss — A/Dlss is the angular size 
corresponding to a comoving distance A at the last scattering 
surface (LSS), whi ch is at a distance D lss — 2(c/.ffof^; 4 ) 
for a fiat universe flVittorio fe Silk 1992; ). Thus 



I ~ k/D LSS ^ 3000fcfi„ 



(10) 



For the F = 0.5 case (corresponding e.g. to f2 m = 0.805, 
SIa = 0.195) the inferred primordial power in Figure |l| is 
suppressed by a factor > 2 for k > 0.1 /iMpc -1 , implying a 
similar decrease in the CMB anisotropy at I > 500. For the 



2.4 Fits to the CMB and LSS data 

We parameterise the 'step' in the primordial power spectrum 
(see Figure [J) as: 



P\k) = 




(ii) 



where k s and k a mark the start and end of the break from a 
H-Z spectrum, with amplitudes A and B. (The values of C 
and a are specified by the other parameters.) In the multiple 
inflation model (Adams et al. 1997b), the actual form of the 
spectrum during the phase transition is difficult to calculate 
since the usual 'slow-roll' conditions are violated. However 
a robust expectation is that 



0.5 < In ( 7r ) < 2, 



(12) 



© 2000 RAS, MNRAS 000, 000-000 



On the APM power spectrum and the CMB anisotropy 5 



because the field undergoing the symmetry-breaking phase 
transition evolves exponentially fast to its minimum. The 
ratio of the amplitudes A/B is determined by the (unknown) 
superpotential couplings of the field undergoing the phase 
transition but is expected to exceed unity (i.e. there is a 
decrease in the power). 

We allow the spectral and cosmological parameters to 
range rather widely as follows: 

• k B from (0.01-0.15) AMpc -1 

• ln(ft e /fc„) from 0.01-4 

• A/B from 0.3-7.2 

• J1a from 0-0.7 (keeping f2 m + £7 a = 1) 

• h from 0.45-0.9 

• Q B h 2 from 0.0014-0.033 

We also consider a bias parameter b, defined as the square 
root of the ratio of the APM and CMB normalisations 
(which is expected to be close to unity, see Appendix A). 

For each combination of the above parameters we con- 
struct the primordial power spectrum and then obtain the 
matter power spectrum by convoluting with the CDM trans- 
fer function calculated by CMBFAST (which differs slightly 
from the analytic approximation in Equation ^). It is also 
used to calculate the expected CMB angular power spec- 
trum. These are then compared, respectively, with the APM 
data ( iGaztanaga fc Baugh 199&j) PI and the decorrelated 
COBE ( Tegmark fc Hamilton 1997 ) and Boomerang data 

( ie Bernadis et al. 200C ) ^1 using appropriate window func- 



tions, to determine x' 2 ( a ) — ^ . [(model(a) — data(i))/cr(j)] 2 . 
We use 19 CMB data points [7 from COBE (excluding the 
anomalously low quadrupole) and 12 from Boomerang] and 
10 APM data points (for 0.01 /iMpc" 1 < k < 0.2/iMpc _1 ). 
We quantify the goodness-of-fit following the procedure of 
Lineweaver & Barbosa (1998). 

In Table |l| we show the values of the Hubble constant, 
baryon density, and amplitude change at the spectral break 
which minimise the \ 2 (f° r a specific choice of Qa) for fits 
to the CMB data alone. The step in the primordial spec- 
trum is seen to be not essential for a good fit, in particular 
even with A/B — 1 (i.e. no step) the x 2 does not worsen 
significantly in most cases. The best overall fit obtains for 
Qa = 0.5, con sistent with the analysi s by the Boomerang 
collaboration ( |dc Bernadis et al. 200C ). Note however that 
the baryon density required is 65% higher than the preferred 
BBN value (|s|) although still within the upper limit (jf]) set 
by interstellar deuterium. 



^ We use the values of Papm(^) given in Table 2 with errors 
estimated from the variance in the 4 disjoint sky regions in the 
catalogue. Ideally we should use the linear spectrum Pu n (k) recov- 
ered for each specific choice of cosmological parameters. However 
the two are very similar since we restri ct o urselves to the quasi- 
linear range k < 0.2/iMpc — 1 (see Fig |Al| in Appendix A), and 
the estimated errors in Papm(^) in this range are large enough 
to encompass all possible linear spectra. 

II There appear to be systematic differences between the 
Boomerang and MAXIMA data as shown in Figure ^. However 
according to Jaffe et al (2000), the two datasets are quite con- 
sistent within their respective calibration uncertainties (20% and 
8%, la in C;) and beam uncertainties (10% and 5%). For definitc- 
ness we use the Boomerang data alone; clearly our results would 
be very similar if we use the MAXIMA data instead. 
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6.6 (14) 


0.4 
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0.031 


> 1.1 


6.3 ( 8) 
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> 1.0 
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0.02 


0.7 
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5.5 ( 6) 
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Table 1. Parameters for best fits to CMB data (13 d.o.f.). 
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0.8 
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6.9 




0.4 
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0.03 


0.8 


1.7-3.1 


6.9 




0.5 
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0.04 


0.6 


1.7-2.8 


7.9 




0.6 


0.75 


0.05 


0.6 


1.5-2.1 


11.9 




0.7 


0.80 


0.01 


0.6 


0.7-0.9 


35.1 





Table 2. Parameters for best fits to CMB data with the BBN 
constraint on the baryon density. 

If we now demand that the models should have the BBN 
baryon density (jE|), the Hubble parameter in the range (Q) 
and a spectral break duration (|l2|), then a step in the pri- 
mordial spectrum is definitely indicated as seen from Ta- 
ble ^. [In particular setting A/B = 1 (no step) always gives 
an unacceptably high value of x 2 > 26.] The angular power 
spectra corresponding to some of these models is shown in 
Figure |§] along with the Boomerang and COBE data. We 
emphasize that a good fit is now obtained even for Oa = 0. 

Next we fit simultaneously to both the APM and CMB 
data allowing for a small difference in the normalisations, 
i.e. a bias b. We have taken into account that the mean red- 
shift of the APM galaxies is z — 0.15 so a corresponding 
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0.60 
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2.4-3.6 


9.0 
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0.023 
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30.4 



Table 3. Parameters for best fits to CMB+APM data (22 d.o.f.). 
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Figure 3. CMB angular power spectra for the best-fit models in 
Table ^| with Qa = (continuous line), fl\ = 0.2 (dashed line), 
Qa = 0.4 (dot-dashed line), and Qa =0.6 (dotted line). The data 
shown are from COBE and Boomerang. All models have the BBN 
baryon density Hb = 0.019/i — 2 . 
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Figure 4. Primordial power spectra for the best-fit models in 
Table ^ with Ca = (continuous line), Ca = 0-2 (dashed line), 
tt A = 0.4 (dot-dashed line), and fl A = 0.6 (dotted line). All 
models have the BBN baryon density f2e = 0.019/i — 2 . 



n A 


h 


k s 




A 
B 




x 2 


0.0 


0.50 


0.07 


1.5 


4.3-6.5 


0.65-0.73 


11.9 


0.2 


0.55 


0.06 


1.6 


3.6-5.6 


0.69-0.78 


10.0 


0.3 
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10.3 


0.4 


0.65 


0.05 


1.5 


3.1-4.6 


0.77-0.87 


13.1 


0.5 


0.70 


0.05 


1.6 


2.6-4.1 


0.81-0.93 


18.8 


0.6 


0.75 


0.05 


2.0 


2.0-2.8 


0.90-0.99 


29.7 



6000 




Table 4. Parameters for best fits to CMB+APM data with the 
BBN constraint on the baryon density. 

correction for the growth factor to z = needs to be ap- 
plied. It is seen from Table ^| that the best overall fit is still 
obtained for Q\ — 0.5, although other values remain quite 
acceptable. In particular the £Ia = model requires a Hub- 
ble parameter at the low end of the allowed range but has 
the advantage of being consistent with the BBN value of the 
baryon density. We emphasise that a break in the spectrum 
is now definitely required since the value of \ 2 exceeds 45 
otherwise. The most likely duration of the break also comes 
out in accord with the theoretical expectation (p^). 

Finally Table ^ shows the result of imposing the BBN 
constraint (^), and the (less important) constraints on the 
Hubble parameter (Q) and the duration (|l2|). We have also 
required that the bias b be within 20% of unity (see Ap- 
pendix A). We see that the data now prefer lower values of 
Qa (and h). The magnitude of the required spectral break 
decreases with increasing Qa (Figure ^J) but one cannot do 
without such a break altogether since otherwise the ex- 
ceeds 53. Figures |E] and ^| show, respectively, the fits to the 
CMB and APM data for these models (using the mean value 
of the break amplitude A/B). In Figure [7J we display the 
goodness-of-fit contours in the fish 2 — h plane for the joint 
fit to the CMB and APM data; the contours are ragged 
because we have not used a fine enough grid of parameter 
values. We see that agreement with the BBN value (Eh of the 



I , i , , , i , I 

200 400 600 800 1000 

multipole 

Figure 5. CMB angular power spectra for the best-fit models 
in Table |j] with Ca = (continuous line), £7a = 0-2 (dashed 
line), Qa = 0.4 (dot-dashed line), and Qa = 0.6 (dotted line). All 
models have the BBN baryon density Ob = 0.019h — 2 . The data 
points are from COBE and Boomerang. 

baryon density can be achieved only for low Qa- In fact a 
critical density matter-dominated universe is favoured with 
the Hubble parameter at its lower limit (so that the age of 
the universe is acceptable at ~ 13 Gyr). 

2.5 Implications for cluster abundances 

In Table ^ we show also the value of the variance o(R) in 
(dark matter) fluctuations (normalised to the CMB), over a 
sphere of size R = 8h~ Mpc: 

1 f°° 

a 2 (R) = -j / W\kR) 5l{k) T 2 {k) k 3 dk, (13) 
■"o Jo 

using a 'top hat' smoothing function, W(kR) = 3[ S ^^P — 
Tfcfl)^ ]' (The range of erg for a given value of Q A cor- 
responds to the range of the step-size A/B.) The tabu- 
lated values are systematically somewhat higher than the 
mean value as = 0.56f^ 0,47 for flat cosmologies inferred 
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Figure 7. Goodness-of-fit contours for the joint fit to the CMB and APM data of a primordial spectrum with a 'step', for different 
choices of Oa- The vertical and horizontal shaded bands indicate the constraints (Q) and (jjj) on h and S1b^ 2 respectively while the crosses 
indicate the best fits from Table |^. Values of f2 A > 0.6 provide an unacceptably poor fit to the data. 



from the observed abundances of X-ray emitting clusters, 
although within the 95 % c.l. upper limit whic h is higher by 
20S^ 21ogf2m per cent ([Viana fc Liddle 1999|). We wish to 



emphasise that the break in the primordial spectrum im- 
plies a much lower as than the corresponding H-Z spectrum 
with the same cosmological parameters. For example the 
COBE-normalised sCDM (H-Z) model with fi m ~ 1 gives 
as — 1.4 (e.g. Stompor, Gorski and Banday 1996) as op- 
posed to as — 0.65 — 0.73 in Table [l|, which is much closer to 
the value of as — 0.56 inferred from the cluster abundance. 
As shown in Figure ^, this is because of the decrease in power 
at small scales relative to a scale-invariant spectrum. Note 
that we cannot constrain well the region k > 0.2/iMpc -1 
(corresponding to R < 8/i~ 1 Mpc), because of uncertain- 
ties in the reconstruction of the primordial spectrum due 
to non-linear biasing and non-linear gravity effects. For sim- 
plicity we have assumed a H-Z spectrum for k > 0.6 h Mpc -1 
(Equation |^) so in principle other (smaller) values of as are 
also allowed. For instance, if there were further breaks in 
the primordial spectrum, as is quite possible, then a better 
agreement could easily be achieved. In order to constrain 
such possibilities we intend to investigate in further wo rk the 
Lyman -q f orest P(k) and its evolution with redshift ( Croft 
et al. 1)99) which provides a probe of such small scales in 
the linear regime. 



In Figure |^ we show the predictions for the redshift evo- 
lution of the cluster abundance in our models according t o 
the Press-Schechter (P-S) theory ([Press fc Schechter 197^), 



compared to the observations as presented in Bahcall & Fan 
(1998). We consider clusters with mass M > 8 x 10 14 M o 
(neglecting the R < R CO m = 1.5 /i -1 Mpc condition). The 
spatially flat COBE-normalised H-Z model with £2a = 0.2, 
Q. m ~ 0.8 (top thick full line) does not match the obser- 
vations, a fact that has been used to rule out this model. 
Even when this model is normalised at z = by lowering its 
amplitude to as — 0.6 (bottom thick full line), it still fails 
to reproduce the observed redshi ft evolution in the cluster 
abundances (Bahcall fc Fan 1998). By contrast, the best-fit 



models in Table |4|do match the redshift evolution, especially 
if we allow for the possibility mentioned above of lowering 
the erg values further to better match the z = data.|^j We 
emphasise that these predictions have not been adjusted in 
any way to fit cluster abundances; all parameters have been 
fixed already by the fit to the APM and CMB data. It is 

** Note that the P-S prediction depends not only on the value 
of as at z = and the corresponding cosmological evolution of 
the variance, but also on the shape of the spectrum over the 
mass scales considered. In our case, these scales correspond to 
R > 10h — 1 Mpc, i.e. near the break in the primordial spectrum. 
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els, the break in the H-Z spectrum discussed here is gener- 
ated by a physical mechanism, viz. supersymmetry breaking 
during inflatio n (Adams et al. 1997b, Ross 1 9 98). As noted 




subseque ntly (Lesgourgues & Polarski 1997 



Kanazawa et 



2000), the resulting damping of the CMB anisotropy 
is also a characteristic of 'double inflation' which occurs 
in toy models wherein inflation is driven in two stag es — 
first by higher- order gravity and then by a scalar field ( | lvof- 



ni.in, Linde & Starobinsky 1985 ; Kofma n & Pogosyan 198S ; 

Starobinsky 199l|), or by two coupled 



Gottlober, Miiller 



scalar fields (Starobinsky 1985 ; Kofman fc Linde 1987 ; Po- 



larski fc Starobinsky 199$ . The model parameters can be 



0.01 



0.10 



k (h/Mpc) 



tuned to generate a break in the primordial spectrum at 
any chosen scale. The advantages of this for fitting LSS 
observations was first emphasised by Silk fc Turner (1987) 
and have been explored thoroughly (|Polarski fc Starobin- 



111 M A 



Figure 

ble | w 

n A = ( 

models have the BBN baryon density Hb 
are from the APM survey. 



sT^frtTnocTcTsTrm^ 
: U.2 pUBSHeT! 11111!), 



sky 1992;; iGottlobcr, Miickct fc Starobinsky 1994 iPolarski 



U {(JollllUUolia 11111!), M A 

4 (dot- dashed line), arid fi A — 0-0 (dotted line). All 



1994; Peter, Polarski fc Starobinsky 1994; Amendola et al 



1995; Polarski fc Starobinski 1995; Semig fc Miiller 1996) 



0.019/1" 



The data 



although these authors did not discuss the implications for 
the CMB. A spectral feature similar to that observed can 
also be generated in a toy mo del where the infla ton evolves 
through a kink in its potential ( [Starobinsky 1992 ), lea ding to 




similar suppression of the secondary acoustic peaks (Atrio 



Barandcla et al. 1997; Lesgourgues, Polarski fc Starobinsky 



1996 ; Lesgourgues, Piunct fc Polarski 1999^ 7 Recently there 



have been efforts to progress beyond toy models and imple- 



ment multiple inflation in supersymmetric models (Sakel- 



lariadou & Tetradis 1998; Lesgourgues 199£; Lesgourgues 



O 
CO 

A 



: r 
o 



200C). Spectral features can also be generated by resonant 



production of particles during inflation (Chung et al. 2000) 



0.4 0.6 
z T red s h if t) 



Figure 8. The number density of rich clusters with M > 8 X 
10 14 A/q as a function of redshift. Points with errorbars corre- 



spond t ) the observations as depicted in Bahcall & Fan (1998). 
The thick mil lines correspond to the P-S prediction lor the spa- 
tially flat, COBE-normalised H-Z model with Q A = 0.2, fl m ~ 0.8 
(top line), and scaled to match the cluster abundance at z = 
(bottom line). The curves in between correspond to the best- fit 
models in Table ^ with £7a = (full line), S7 A = 0.2 (dashed line), 
Qa = 0.4 (dot-dashed line), and Q\ = 0.6 (dotted line). 



interesting to note that even the J1a = model can now 
reproduce the slow growth of the cluster abundance with 
redshift due to the break in the primordial spectrum. 



Silk & Gawiser (2000) have noted independently that 
it is difficult for a ACDM model to fit LSS and CMB ob- 
servations simultaneously without considering departures 
from scale-invariance; they advocate adding a 'bump' at 
k ~ 0.06 /iMpc -1 to a primordial H-Z spectrum.^] Einasto 
et al. (1999) have also concluded from an exhaustive analy- 
sis of LSS data that a scale-free primordial spectrum is ex- 
cluded. In particular it has been noted that the power spec- 
trum of clusters is better fitted if the pri mordial spectrum 
has a step- like feature at k ~ 0.1 ftMpc -1 (Gramann & Suh 



honenko 1999). The implications of broken scale-invariance 



for accurate determination of cosmological parameters us- 
i ng forthcoming LSS and CMB data have been in vestigated 
( |Wang, Spergel fc Strauss 1999| [Hannestad 2000| ). 



It is clear that the notion of a scale imprinted on the 
primordial density perturbation is favoured by a number 
of observational considerations as well as having a natural 
physical interpretation in the context of realistic inflationary 
models. In this paper we have demonstrated that the recent 
small angular-scale CMB anisotropy data favour a deviation 
from scale-invariance in the range k ~ (0.06 — 0.6)/iMpc -1 
independently of the similar previous indication from stud- 



3 DISCUSSION 

Several authors have in the past noted the possibility of 
generating spectral features in the H-Z density perturba- 
t ion by fine tuning the parameters of the inflaton potential 



(|Kofm£jn~fc Linde 19871; kalopck, Bond fc Bardeen 198^ 



HodgesTFal. 1990| ; [Hodges fc Blumcnthal 1990| ). We woukf 



like to emphasise that in contrast to such "designer" mod- 



tt Griffiths, Silk &: Zaroubi (2000) have very recently shown that 
this can also fit the recent CMB and LSS data; however they 
choose to discount the APM data in favour of the ("le ss nre- 
cisei power spectrum inferred from the PSCz survey (Hamil- 
ton & Tegmark 2000) which does not show a 'shoulder' at 



k ~ 0.1 /iMpc T . We would point out that our models in Ta- 
ble ^ also fit the PSCz power spectrum. 
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ies of LSS, in particular the APM survey. This is reassur- 
ing since the latter rests on the assumption (for which ar- 
guments are given in Appendix A) that there is no scale- 
dependent bias between APM galaxies and dark matter. 

Allowing the primordial spectrum to depart from scale- 
invariance has important consequences for the deductions 
that can be made about cosmological parameters from LSS 
and CMB data. The standard interpretation of LSS data 
(e.g. APM P(k) and erg from cluster abundances) favours a 
low r ~ Q m h ~ 0.2 universe assuming a COBE-normalised 
H-Z spectrum. However as we have shown these observa- 
tions can be satisfactorily accounted for with much larger 
values of Q m if we allow for a break in the primordial spec- 
trum. Our best fit models in Table ^, which obey the BBN 
constraint (||) on the baryon density, prefer higher values of 
Qm (and lower values of Qa)- The required amplitude of the 
break decreases with increasing Qa but one cannot do with- 
out such a break altogether. In particular the value Qa ~ -7 



favoured by the SN la data (Perlmutter et al. 1999; Riess et 
al. 199^) is not permitted, while it is quite possible to have 
a universe with zero cosmological constant and Q m = 1. 
Of course the high observed fraction of baryons in clusters, 
combined with the BBN value of the baryon densi ty 
implies an independent limit on the matter density (White 
et al. 1993 ); a recent analysis quotes Q m < 0.31 h~ llz at 
95% c.l. (|Mohr, Haiman fc Holder 200C|). There are many 



assumptions made in such analyses, viz. that clusters are 
spherical and in hydrostatic equilibrium with an isothermal 
temperature profile, that the measured baryon fraction re- 
flects the universal value, that there is no preheating before 
or energy injection after cluster formation, etc. It is impor- 
tant to reassess these issues in the light of our increasing 
understanding of galaxy and cluster formation to quantify 
the systematic uncertainties better. 

The hypothesis of a primordial density perturbation 
with broken scale-invariance, although radical, is eminently 
falsifiable. The ongoing 2DF and SDSS redshift surveys can 
confirm or rule out such a feature in the power spectrum of 
galaxy clustering, while the forthcoming MAP and Planck 
missions will determine whether all the secondary acoustic 
peaks in the CMB angular spectrum are indeed suppressed 
as expected. By contrast the alternative hypothesis of a 
baryon density ~ 65% higher than the BBN value predicts 
a boosted third peak (e.g. Tegmark & Zaldariagga 2000), so 
this would be a definitive test distinguishing between the two 
possibilities. We note that in either case there are important 
implications for the physics of the early universe. To increase 
the baryon-to-photon ratio after BBN by ~ 65% requires a 
decrease in the comoving entropy; the only mechanism sug- 
gested which can achieve this proposes that the photons 
in our universe be come cooled through their coupling to a 
'shadow' universe (Bartlett & Hall 1991). However this pos- 



sibility, recently invoked by Kaplinghat and Turner (2000), 
was shown to be severely constrained, if not ruled out al- 
together, from considerations of the concommitant spectral 
distortion in the CMB together with the bound fro m Su- 
pernova 1987a on new weakly interacting particles ( Birke] 



Sarkir 1997). By contrast broken scale-invariance has a 



natural explanation in a phase transition occuring during 
inflation as expected in supersymmetric theories. If estab- 
lished this would provide the first direct connection between 
astronomical data and physics at very high energies. 
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APPENDIX A: THE EFFECTS OF BIASING 

Independently of how or where galaxies form they must 
eventually fall into the dominant gravitational wells and 
thus trace out the underlying mass distribution (see Peebles 
1980, Fry 1996). However, in principle, the galaxy distribu- 
tion may be 'biased' with respect to the underlying mass 
fluctuations (e.g. Bardeen et al. 1986) . There are several ob- 
servations indicating that this effect is in fact small. The 
theoretical predictions for the first few connected moments, 



based on this hypothesis ( 


Juszkiewicz, Bouchet & Colombi 


1993; 


Bernardeau 1994a; Bcrnardcau 19994b|) arc in good 


agreement with the APM data (Gaztanaga 1994; 


Gaztanaga 


k Frieman 1994 




Frieman & Gaztanaga 1999). The current 



precision of this higher order correlation test is 20% and ex- 
pected to improve with further data. The APM angular cor- 
relation function £ g (r) shows an inflection point near £ g ~ 1 
very similar to the long awaited ' shoulder' in £ ~ 1 , gen- 
erated by gravitational dynamics flGott fc Rees 1975[ ). The 



agreement between the two characteristic scales can be used 
to constrain the linear bias ing factor for the APM catalogu e 
to be within 20% of unity ( Gaztanaga fc Juszkiewicz 200C ). 
We summarise these arguments ( Gaztanaga & Baugh 199c ) 



below for completeness since the assumption of small or no 
bias is crucial for the present study. 

Let us assume that the (smoothed) galaxy fluctuations 
S g are related to those in the mass S m by a local transforma- 
tion, 8 g (x) — F[8 m (x)], which can be expanded as a Taylor 
F — bi8 m + b-zS'm + . . .. Then the 2-point function 
■ (8g(x)S g (x + r)) on the scale r will be: 



series 



(r) 



CIM = 6?£ 2 m (r) + &i& 2 (Mx)Mz + r) 2 ) 
+bib2(dli(x)5 m (x + r)} + 



(Al) 



where all further terms are of order 4 or higher in S m and 
correspond therefore to either higher order correlations, £j 
with J > 2, or higher powers in £2- If 8 m is Gaussian or 
hierarchical (as is the case for evolution through gravity), 
the higher order correlations £j are at most of order ■ 
This means that at large scales where £2 < 1, the first term is 
dominant so that only the amplitude of the 2-point statistics, 
but not its shape, may be altered by biasing. This has been 
confirmed in N-body simulations and galaxy bias ing models 
( |Mo, Jing fc White 1997| ftcoccimarro Tet al 2000| ). 

Thus from the above arguments, the small variance on 
large scales (R > 8/i _1 Mpc) means that it is reasonable 
to assume that the shape of P(k) for galaxies at scales 
k < 0.1/iMpc _1 coincides with the shape of the underly- 
ing linear matter power spectrum. This argument, which 
is based simply on the smallness of the variance and the 
assumption of hierarchical clustering, can also be applied 
to gravity, since the leading contribution to the correlation 
functions in perturbation theory is indeed given by a lo- 
cal transformation (see Fosalba & Gaztanaga 1997). This 
is clearly illustrated in Figure 11 of Gaztanaga and Baugh 
(1998). By comparing the linear and non- linear shape of 
P(k), one can see that it has not been changed significantly 
through gravitational evolution on scales where the rms fluc- 
tuations are small, i.e. for k < 0.1 ftMpc" 1 . 

How important are non-linearities for 8 > 1, i.e. for 
1? 



Figure 



Al 



shows the ratio of the non- 



k > 0.1/iMpc" 
linear to linear power spectrum P(k), normalised to be 
unity at small k in order to emphasize the non-linear ef- 




0.05 



0.1 

k (h/Mpc) 



0.5 



Figure Al. The ratio of the non-linear to linear power spectrum 
illustrating non-linear effects on P(k). Long and short-dashed 
lines show non-linear effects in the mass while the other lines 
show non-linear effects in the final galaxy distribution for various 
halo models (from Scoccimarro et al 2000). 



fects. The long-dashed line corresponds to the APM-like 
P(k) for Q. = 1 dynamics (from N-body simulations in 
Baugh & Gaztanaga 1996). The short-dashed line corre- 
sponds to a Q = 0.3, A = 0.7 ACDM model (using the 
non- linear mass fitting formulae). The dotted, continuous 
and dot-dashed lines correspond to different non-linear mod- 
els for galaxy formation based on halo models (with linear 
bias b < 1.5 for trg ~ 1) from Figure 6 in Scoccimarro et 
al (2000). In all cases non-linear effects become significant 
(> 10%) only at k > 0.2/iMpc^ 1 and are of order unity by 
k ~ 1/iMpc -1 . Non- linear effects seem smaller for galaxies 
than for the dark matter (in agreement with the observed Sj 
in Baugh & Gaztanaga 1996). One could imagine construct- 
ing galaxy formation models with stronger non-linear ef- 
fects, but it will then be difficult to account for the observed 
galaxy higher-order correlations. As pointed out in Scocci- 
marro et al (2000), galaxy clustering in the framework of the 
halo models is strongly constrained by the measured galaxy 
skewness. Thus, given current observational constraints and 
present understanding of galaxy formation one can conclude 
that non-linear effects cannot have altered the linear P(k) 
by a factor exceeding ~ 2, up to k ~ 1 ZiMpc -1 . 
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